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INTRODUCTION 


It  is  known  that  when  projectiles  are  fired  from  gun  tubes,  the  lateral 
forces  of  constraint  may  become  so  large  as  to  cause  engraving  of  the  projectile 
body  and  excessive  wear  of  the  rifling  of  the  tube  -  usually  at  the  muzzle  where 
the  projectile  velocity  is  highest,  thus  giving  rise  to  the  term  muzzle  wear. 
Normally,  a  projectile  will,  at  some  time  during  its  travel  down  the  bore  of  a 
cannon  tube,  come  to  bear  against  a  particular  rifling  land  and  follow  that  land 
out  to  the  muzzle;  and  when  the  same  land  is  borne  against  with  great  enough 
force  by  a  large  enough  number  of  projectiles,  the  result  may  be  what  is  some¬ 
times  referred  to  as  spiral  wear.  In  any  case,  the  result  is  undesirable  and  we 
wish  to  understand  it  so  that  it  can  be  avoided. 

All  projectiles  do  not  necessarily  ride  on  the  same  land  for  the  full 
length  of  travel  down  the  bore,  and  when  the  balance  or  alignment  of  a  projec¬ 
tile  is  such  or  becomes  such  that  the  shift  from  one  land  to  another  is  violent 
enough  to  cause  rebound,  then  the  phenomenon  called  bal lotinq  may  have  begun. 
Since  it  is  widely  held  that  projectile  balloting  can  have  serious  consequences, 
it  has  been  the  subject  of  much  study,  initially  by  Reno  (ref  1),  and  Thomas 
(ref  2),  and  more  recently  by  Perdreauville  (ref  3),  Chu  and  Soechting  (ref  4), 
Walker  (ref  5),  and  Soifer  and  Becker  (ref  6).  Although  it  is  not  t.ie  purpose 
of  this  study  to  go  into  the  subject  of  balloting,  there  may  be  some  application 
of  the  results  to  the  problem  of  determining  the  necessary  conditions  of  projec¬ 
tile  balance  and  alignment  and  tube  motion  and  curvature  for  balloting  to  occur. 

The  problem  that  is  addressed  here,  however,  is  rather  that  of  determining 
the  combined  effect  of  projectile  balance  and  alignment  and  tube  motion  and  cur¬ 
vature  on  the  bearing  (constraint)  forces  between  the  projectile  and  the  tube 

References  are  listed  at  the  end  of  this  report. 
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when  stable  projectile  motion  takes  place.  Darpas  (ref  7)  used  some  simplifying 
approximations  in  order  to  isolate  the  various  effects  and  estimate  their  rela¬ 
tive  maximum  magnitudes.  His  conclusion  was  that  the  dynamic  unbalance  of  the 
projectile  would  most  likely  have  the  greatest  effect.  Dynamic  unbalance  would 
occur  when  a  projectile  becomes  cocked  in  the  bore,  which  as  discussed  by  Gay 
(ref  8)  and  Montgomery  (ref  9),  is  the  more  usual  case.  The  following  six- 
degree  of  freedom  treatment  of  the  problem  is  intended  to  establisn  the  magni¬ 
tudes  of  the  various  effects  of  specified  projectile  and  tube  conditions. 

ASSUMPTIONS  OF  THE  MODEL 

The  projectile  is  constrained  by  the  inner  wall  of  the  tube  and  the  rifling 
t1  move  in  the  direction  of  and  rotate  about  the  axis  of  the  tube.  In  general, 
the  principal  axis  of  inertia  of  the  projectile  is  skewed  with  respect  to  the 
spin  axis.  The  gas  pressure  driving  force  acts  along  the  tube  axis,  the  land 
driving  force  acts  around  the  circumference  of  the  bore  and  perpendicular  to  the 
rifling,  ana  the  lateral  constraint  forces  act  at  the  point  of  contact  of  the 
rotating  band  and  bourrelet  with  the  bore  and  in  a  direction  perpendicular  to 
the  bore  axis. 

In  addition,  there  are  three  frictional  resistance  forces  normal  to  the 
Isnd  driving  2nd  iB’tsrsI  cons’trsin't  forces  2nd  opposite  to  ths  direction  of  pro- 
jectile  motion  at  the  point  of  contact.  The  band  engraving  resistance  force 
which  acts  around  the  circumference  of  the  bore  and  opposite  to  the  direction  of 
projectile  motion  is  assumed  to  be  known,  as  is  the  gas  pressure,  axial  veloc¬ 
ities  of  projectile  and  tube,  mass,  center  of  gravity  (c.g.)  location,  axial  and 
transverse  moments  of  inertia  of  the  projectile,  and  coefficients  of  friction 
between  the  rotating  band  and  the  tube,  and  between  the  projectile  body  ana  the 
tube. 
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Three  force  equations  written  with  reference  to  the  inertial  coordinate 
system  and  three  moment  equations  written  with  reference  to  the  body-fixed  coor¬ 
dinate  system  constitute  the  six  equations  used  to  solve  for  the  six  unknowns  of 
the  problem.  The  six  unknowns  are  the  two  components  of  each  of  the  two  lateral 
forces  of  constraint  (four  unknowns),  the  driving  torque  around  the  spin  axis 
(one  unknown),  and  the  acceleration  in  the  direction  of  the  spin  axis  (one 
unknown) . 


rB  ■  S{3  +  Pb 


(5) 


Texts  on  dynamics  (ref  10)  give  the  formula 

••  ••  •  • 

pe  3  pb  +  ub/B  *  Pb  *  2u)b/e  x  Pb  *  <*>b/B  *  <*>b/B  *  Pb  <6) 

for  rotating  coordinate  systems  where  the  subscript  b/B  indicates  rotation  of 
the  b  coordinate  system  with  respect  to  the  B  coordinate  system.  If  em  is  the 
magnitude  of  the  eccentricity  of  the  projectile  c.g.  with  respect  to  the  spin 


axis,  then  pb,  pb,  and  pb  can  be  expressed  as 

Pb  =  em  (cos9j  +  sinQk)  , 


and 


pb  ■  emfl  (-sinfij  cos9k)  , 


Pb  =  en[-(0  sinfl  +  0*  cos0)j  +  (0  cose  -  0*  sin0)k] 
where  0  is  as  shown  in  Fiaure  2. 


(7) 

(7a) 

(7b) 


Figure  2.  Position  vector  for  the  projectile  c.g.  in  the  moving 
coordinate  system. 

Since  the  spin  axis  is  taken  in  general  to  be  the  tangent  to  a  curve  moving 
in  space,  it  will  be  convenient  to  refer  to  the  Frenet-Serret  Formulas  of  dif¬ 
ferential  geometry.  These  are 


dB  - 
—  s  -tN 
as 


dj 

ds 


s  kN  ,  and 


5 


dN 

ds 


tB  -  kT 


(8a, b,c ) 


where 


rip  1 

B  =  T  x  N.  ,  T=^,and 

as  P 

The  velocity  of  a  point  moving  along  a  curve  in  space  is  expressed  as 


vit)  =  ^  =  $2.  dR 
y  '  dt  dt  ds 


(9) 


Making  use  of  the  Frenet-Serret  formulas  allows  the  velocity  to  be  written  as 


;<t>  =  t  , 


(10) 


and  the  acceleration  to  be  written  as 


(11) 


Now  since  Sg  is  the  sum  of  the  acceleration  of  the  origin  of  the  b  coor¬ 
dinate  system  relative  to  its  path  and  the  acceleration  of  the  point  on  the  path 
coincident  to  the  b  coordinate  system  origin,  we  can  write 

SB  *  Qb  +  T  +  (~)  <CN  .  (12) 

Also  note  that  w^/g  is  the  sum  of  the  angular  velocities  due  to  rotation  of  the 

spin  axis  at  the  point  coincident  with  the  origin  of  the  b  coordinate  system  and 

motion  of  the  coordinate  system  origin  relative  to  its  path.  Therefore,  if  we 

designate  by  the  subscript  b',  a  coordinate  system  wh-ch  is  coincident  to  the  b 

coordinate  system  but  fixed  to  the  spin  axis,  we  can  write 

■+  -*■ 

“b/B  *  «b'/B  +  “b/b'  -  (13) 

Recalling  that  has  no  component  in  the  direction  of  the  spin  axis,  it 

follows  that 

-  ds  "  dT  ds  ■* 

"b/b*  =  dtTx3i  =  3tTx!cN  (14) 


6 


I 


ds 


so  that 

"b/8  =  wb</g.=  +  -j=  T  x  kN 
Differentiating  with  respect  to  time  gives 


-  ds  P 

wb'/e  +  at  KB 


(15) 


d  .ds  r, 
"b/B  -  «b'/B  +  at  [dt  *B] 


(16) 


which,  after  some  manipulation,  results  in 


.  d*s  die  .ds,*  r  .ds,p 
“b/B  ■  “b'/B  +  t*  at*  +  as  (dt)  38  “  *T(dt)N 


(17) 


COORDINATE  TRANSFORMATION 

^  **♦ 

First  of  all,  by  observing  that  T  ■  i,  and  N  =  Nyj  +  Nzk,  it  follows  that 

-»  *♦  *+  -*•  ~+ 

B  a  TxN  a  Nyk  -  Nzj.  Also, 


and 


Qb  3  (QB)xi  +  (QB)yJ  +  (QB)zk  > 

-> 

1 /B  =  (ub'/B)yJ  +  (Wb'/B)zk  < 


(18a) 

(18b) 


«b'/B  *  (“b'/BJyJ  +  (wb«/B)zk 


(^c) 


Therefore, 


"*  ”*  d*s  ds  *  ”** 

SB  *  [(Qb)x  ♦  3t*^i  ♦  C (Qb )y  +  (3^)  KNy^J 


♦  hqb)z  ♦  (af)*<Nz]k 


“b/B  *  t(«b’/B)y  -  ai  kNz]j'  +  t(«b'/B)z  +  a*  KNy]*  » 


(19a) 

(19b) 


and 


r  i  *  ,  .  d*s  die  ,ds  *...  .ds  *  ip 

«b/B  -  {(«b'/B)y  -  t*  3ti  +  di  (dt)  3Nz  "  KT{dt>  Ny)J 


« ,  .  d*s  die  ,ds,*1u  ,ds  *  ip 

+  U^'/bJz  +  tK  at*  +  3s  (3t)  ]Ny  ”  <t(je)  Nz)k 


(19c) 


The  total  angular  velocity  of  the  projectile  is  the  vector  sum  of  the  pro¬ 
jectile  spin  velocity  and  the  angular  velocity  of  the  b  coordinate  system,  or 

^  ^ 

n  =  w  +  wb/B  .  (20) 


i 
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In  the  b  coordinate  system,  the  components  of  C  are  then 


flx  *  w  , 

Qy  =  (<*»b'/B}y  +  kN2  •  and 

oz  =  («bve)z +  if  kNv  * 

Similarly,  the  angular  acceleration  and  its  components  are  given  by 

t  •  t 

Q  -  «  ♦  «b/B  » 

•  * 
ax  *  w  , 


and 


°y 


Q-, 


»  r  d*s  die  .ds . 

(Wb'/B)y  “  [K  Jt*  +  ds  *dt*  ^Nz 


.**  .  .  d*s  die  ,ds 

(<^b'/B)z  +  dF  +  ds  (dt}  ]Ny 


ds  * 

KT<dt}  Ny  ' 


KT(5f)aNz 


(21a) 

(21b) 

(21c) 

(22) 

(23a) 

(23b) 

(23c) 


It  remains  now  to  express  the  projectile  angular  velocity  and  acceleration 
in  terms  of  their  components  in  the  body-fixed  coordinate  system.  For  the  angu¬ 
lar  velocity 

Qj  ■  cosfiixQx  +  cosSjyQy  +  cos6jzQz  ,  (24a) 

fi2  =  cos62xQx  +  cosdgyOy  +  cosd2znz  ,  (24b) 

and 

Q3  «  cos63Xflx  +  cos53yt2y  +  oosfi3Zftz  .  (24c) 

Remembering  that  5jx,  62x,  and  $3x  are  assumed  to  remain  constant,  we  get  by 
differentiation 

•  *  d  * 

Qj  *  C0S6jXGx  +  (COSdly)fiy  +  COSdlyfly 


+  <cos6lz>°z  +  cos6lzfiz 


(25a) 


Q2  =  cosd2xnx  +  (cos52y)Qy  +  cos62yfly 
+  ^  (cos d2z)nz  +  cosd2zfiz  . 


(25b) 


-»'•  r  ,  V 


and 


Q3  *  cos«3x0x  +  (cos«3y)Oy  +  cos63yfly 

+  ~  (cos53z)n2  +  cosa3znz  .  (25c) 

Of  the  nine  direction  cosines  appearing  in  the  above  equations,  only  two 
are  actually  independent.  The  other  seven  are  expressible  in  terms  of  these  two 
and  are  therefore  eliminated  as  physical  parameters.  Because  the  conversion  is 
lengthy,  it  is  presented  separately  in  the  Appendix,  and  only  the  results  are 
shown  here. 

W5 

cos6iy  =  sin6jx  cos  —  ,  (26a) 

cos6jz  *  sin6lx  sin  ~j  ,  (26b) 


and 


COS62y  * 


COS&2Z  = 


cos63y  = 


COSSjj,  a 


ns  ns 

-  cos6jx  COS62X  COS  Hi  "  s"'r*  i~N  cos63x 

sinfilx 

ns  ns 

-  cos6jx  cos52x  s’’n  rN  +  cos  rN  cqs63x 

sin6ix 

ns  ns 

-  cos£i..  me  TTi  +  cSn  “  rnQ.W. 

aa - va - i  re  in  — -  ~  4A 

Sin5lx 

ns  ns 

-  cosfijx  cosfi3x  sin  -  cos  JT]5  cos62x 

sin^jx 


(260 


(26d) 


(26e) 


(26f ) 


cos*6jx  +  cos*d2x  +  cos*63x  a  1  .  (27) 

Performing  the  differentiation  indicated  on  the  direction  cosines  in  the 
formulas  for  the  angular  acceleration  components  reveals  the  following: 
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“t(C035ly)  *  -  C0S6lz  (  — )  , 
ft  (cos <5U)  =  cos6ly  (J|)  , 
ft  (cosa2y)  =  -  cos62z  (Jg)  . 

d  .  e  .  *  .7TS. 

—  (cos 622)  »  cosfi2y  (pg)  . 

d  .  _  .  .  .7TS. 

dt  ^cosfi3y)  a  "  cos«3z  (pg)  , 

and 

ft  (cosfi3z)  =  cosfi3y  (pg)  , 

Therefore, 

i  •  ® 

*  COS«lXflX  +  COSdly  [Qy  +  (pg)QZ] 
ire 

+  cosdl2[a2  -  {';^)ny3  , 

•  *  *  ns 

Q2  *  cos52xqx  +  cos62y[fly  +  (j^j)  flz3 
^  • 

♦  cosfi22[nz  -  (pg)Qy]  » 

and 

'  '  '  7T5 

Q3  =  cosfi3xnx  +  cos63y [Qy  +  (pg)nz3 


+  cos63z[£2z 


Sly]  . 


FORCES  ON  THE  PROJECTILE 

Three  of  the  six  equations  reeded  for  solution  of  the  problem  are 
from  the  equation  of  motion 


n 

2  ^  «  mP^B  • 

i*l 


(28a) 

(28b) 

(28c) 

(28d) 

(2Re) 

(28f) 

(29a) 

(29b) 

(29c) 

avai lable 

(30) 
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Of  the  six  forces  acting  on  the  projectile,  one  is  the  known  gas  pressure 
driving  force  and  two  are  the  unknown  lateral  constraint  forces  of  the  problem. 
The  three  remaining  forces  are  frictional  resistance  forces  which  arise  from  the 
band  pressure  force,  spinning  torque,  and  constraint  force  applied  between  the 
projectile  and  tube  in  relative  motion. 


Figure  3.  Forces  in  the  developed  circumferential  (x-9)  plane. 

Figure  3  is  a  development  of  the  rifling  curve  onto  the  plane  defined  by 
the  x  and  0  axes.  The  band  pressure  drag  force  is  denoted  by  D,  the  normal  land 
driving  force  is  denoted  by  Fn,  and  the  angle  between  the  tangent  to  the  rifling 
curve  and  the  x-axis  is  denoted  by  If  we  assume  that  D  is  proportional  to 
the  width  of  the  rotating  band  in  the  direction  of  its  movement  relative  to  the 
tube,  then  we  can  write 

D  =  B/cosi//  .  (31) 

Summing  forces  in  the  x  direction  gives 

Fx  *  "  Fn  sinM?  -  (MBFn+D)cos4>  ,  (32) 

and  since  the  sura  of  the  forces  in  the  6  direction  are  equal  to  the  spinning 
torque  divided  by  the  radius  of  the  bore,  we  have 

Mx 

~  •  Fn  cosip  -  (nBFn+D)sirnjf  (33) 

11 
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Fn  can  be  eliminated  between  Eq.  (32)  and  Eq.  (33)  by  rearranging  and  dividing 


I 


one  by  the  other 


-  Fx  -  0  cos*  sin*  +  mb  cos* 
"5  +  0  sin*  *  cos*  -  MB  sin^ 


Dividing  the  numerator  and  denominator  of  the  right  side  of  the  above  equation 


by  cos  *  gives 


-  Fx  -  B  ^  f 1 (s)  +  mb 
^  ♦  Bf'(s)  "  1  *  ^f'(S) 


For  constant  pitch  rifling  f'(s)  *  r  so  that 


F  .  *  .  Brr  ”  -  MBN  _ 
x  *r  N  *'mb*  _  N* 


Lateral  Force  Friction  Components 


The  lateral  force,  F~  at  the  band  and  the  lateral  force,  Fp,  st  ths 


bourrelet  have  frictional  forces  associated  with  them  which  have  magnitude  mbFr 


and  MsFp,  respectively.  These  frictional  forces  act  along  the  bore  surface 
opposite  in  direction  to  the  motion  of  the  projectile  and  therefore  have  axial 


and  tangential  components  as  shown  in  Figure  4. 


Figure  4.  Components  of  the  lateral  force  associated  friction  forces. 


I 


ft 


ft 


r»s 


I 


i 


Again,  for  constant  pitch  rifling  tar^;  »  f'(s) 


2  ,  so  we  have 
N 


and 


The  tangential  component  can  now  be  broken  up  into  y  and  z  components  as 


Figure  5.  Lateral  force  components  in  the  moving  coordinate  system 


From  the  similarity  of  the  force  triangles  in  Figure  5,  the  following  is 
derived: 


HBFR(n/N)  ( FR )  z  MB<FR)z(n/N) 

■  —  -  -  -  3  —  ~  —  ■  '■  3  —  — — 

/(tt/N)*  ~Tj  /(n/N)*  +  1  Fr  /(tt/N)»  +  1 


(39a) 


tiBFp,(n/N) 
/(n/N)*  +  1 


MBFR(>V/N)  (Fr)v  _  MB(FR)yQT/N) 


✓(n/N)*  +  1  Fr  /(n/N)2  +  1 


(39b) 


-1  z 


and  similarly  for  the  force,  Fp. 


MOMENTS  ABOUT  THE  BODY-FIXED  AXES 


In  addition  to  the  three  equations  derived  from  the  equation  of  motion,  we 
require  three  more  equations  to  give  us  the  six  necessary  for  a  solution.  These 
additional  equations  are  derived  from  Euler's  equations  which  are  written  with 
respect  to  the  body-fixed  axes  as  follows: 

Mj  »  QlIl  -  C2Q3d2-l3)  »  (40a) 

M2  *  02*2  ~  SI3Q1  ( I3" 1 1 )  t  (40b) 


and 


M3  “  Q3I3  -  QjQ2(Il"l2)  •  (40c) 

Note  that  the  quantities  1^  and  I2  are  the  transverse  moments  of  inertia  of  the 
projectile.  Now  since  our  concern  is  with  repeatable  conditions  that  occur  with 
regularity,  there  is  little  point  in  differentiating  between  the  two  transverse 
moments  of  inertia  of  projectiles  which  are  usually  loaded  with  random  angular 
orientetion  and  are  in  any  case  nominally  equal.  Therefore,  we  shall  say  lj  = 

I2  *  It  and  redesignate  the  axial  moment  of  inertia,  I3,  as  1^.  The  Euler 
equations  then  become 
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Mj  -  Qjlj  -  Q2Q3(IT-Ia)  . 


(41a) 


m2  ■  Q2jT  ~  Oa^ldA-1!)  ' 


(41b) 


m3  *  q3*A 


(41c) 


The  moments  about  these  axes  are  given  by  the  triple  product 


Mn  ■  M  •  un  =  (rxF)  •  un  , 


where  r  is  the  position  vector  from  the  origin  of  the  body-fixed  axes  to  the 


point  of  application  of  the  force,  F,  and  the  unit  vectors,  un,  for  the  body- 


fixed  axes  a^e  expressed  in  terms  of  the  direction  cosines  as  follows: 


un  *  cosfinx  i  +  cos6nv  j  +  cos6n2  k  (n=l,2,3) 


For  the  axial  forces,  the  unit  force  vector  is  i  so  that 


_  _  ,  .  ns  7  rs 

r  *  -  pb  =  -  em(sm  k  +  cos  -r.  j)  , 


r  x  i  -  -  •M(sin  K  j  -  cos  ■§  k) 


(rxi) 


(46a) 


(rxi)  •  u2 


cos  6-t 


(46b) 


*•  **  -*  cos 

i,u;i  .  ...  _  _ _ t* 

w~,/  “J*  “  sirTiTv  * 


i  j  \ 

V  «*uc ; 


The  unit  moment  vector  for  the  rotating  band  torque  is  i,  so  that  for  right-hand 


twist  we  have 


*  Uj  •»  cosfiix  , 


(47a) 


i  •  u2  ®  cos52x  > 


(47b) 


i  •  u3  *  cosd3x 


(47C) 


I 


§ 


IS 


10 


1*1* 


K* 


s 


Ki 


y» 


|*K 


l*f* 


la  j.  jk«n  w»intkHjnLA  AjlILa  hat  j  ;vlT^J"a;vcvcv\;  W  A  rfVaV  KV  a-TiJT 


Tha  position  vectors  of  the  points  of  application  of  the  lateral  constraint 
forces  and  their  associated  forces  are  shown  in  Figure  6.  The  distances  from 
the  c.g.  of  the  projectile  along  the  x-axis  to  the  rotating  band  and  the 
bourrelet  are  given  t,  b  and  a,  respectively. 


Figure  6.  Application  points  of  the  lateral  constraint  forces  and  their 
associated  friction  forces. 

From  Figure  6  it  can  be  seen  that 


**F  *  "  p~  [Wyj  +  Wz*]  , 

(48a) 

and 

rR  -  -  j='  C(FR)y3  ♦  (FR)2k]  . 

(48b) 

The  position  vectors  are  then  written  as 


T  tts  “f  *•  -  7TS  •  rp  -  - 

ai  -  [em  cos  +  jrp  (Fp)y]j  -  [em  sin  ^  ^  (FF)2]k  , 

and 

-  bi  -  [em  cos  +  ^  (;R)y]J  -  [6|n  sin  ♦  jjjj  (FR)2jk  . 

Taking  the  vector  product  of  the  position  vectors  and  the  respective  force  vec¬ 
tors  yields  the  moment  vectors 


em C ( Ff )y  sin  p^  -  < Fp )2  cos  pjj  ]i  -  a(FF)2j  +  a(FF)yk  , 
and 

effl[<FR)ysin  2*  -  (fr)z  cos  25]i  +  b(FR)zj  -  b(FR)yk  . 


For  the  sake  of  convenience,  make  the  following  abbreviations: 

Mb 


and 


*»'B 


V{k/h)*  +  1 
Us 


/■{JT/N)*  +  1 

Using  the  above  notation,  the  frictional  forces  associated  with  the  lateral 


forces  are 

M'sC-Fpi  -  <FF)z(J)j  ♦  (FF)y(j[)io  » 
and 

M'Bt'Ffti  "  (FR)z(R)j  +  (FR)y(S)k]  • 

And  the  vector  products  of  the  respective  position  vectors  and  the  above  forces 
are 

-  M's(J)[em(FF)y  cos  p5  +  em(FP)z  sin  p5  +  rFF]i 
-  4’s[o(j-)(FF)y  -  emFp  sin  25  .  r{FF)z]j 

“  H*S[a(jj)(FF)2  ♦  emFP  cos  25  «.  r(FF)y]k  , 

and 

"  cos  Ffi  +  em ( Fr ) z  sin  p5  +  rFp]i 

rf  trc  -*■  •* 

+  M'BCb(fj)  ( Fr )y  +  emFR  sin  p^  +  r(FR)2]j 

+  M'8tb(2)(FR)z  -  eraFR  cos  2®  *  r(FR)y]k  . 

The  moment  of  the  lateral  constraint  forces  and  the  associated  friction 


forces  about  the  origin  of  the  body-fixed  axis  coordinate  system  is  then 
expressed  as: 
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KtWy  sin  PR  "  <FF>z  cos  PR  +  (^R)y  sin  Ffi  '  (fR)z  cos  FR1 

*  s(r>  tem{Fp)y  COS  pjq  +  em(Fp)z  sin  p^  +  rFF] 

“  M'B(R)[em(FR)y  cos  p-q  +  em(FR)z  sin  p§  +  rFR]}i 

+  {-  «(ff)z  +  fe(FR)z  ~  M's[a(R)  ( FF )y  "  *mFF  s-in  PR  '  r<Fp)z] 

*  M'B[b(^)(FR)y  ♦  emFR  sin  p^  +  r(FR)z]}j 
+  {»(FF)y  -  b(FR)y  -  U's[8(rMff)z  +  emFF  cos  pR  +  r < ff )y] 

+  M‘B[b(J)(FR)2  -  emFR  cos  jig  -  r(FR)y]}k  . 

The  components  of  this  moment  vector  in  the  direction  of  the  body-fixed 
coordinate  axes  are  determined  by  taking  the  scalar  products  with  the  unit 
vectors  of  the  body-fixed  axes  cos3nxi  cos^nyj  +  cos$nzk  (n  *  1,2,3)  and 
adding  to  the  axial  force  moment  and  the  spinning-up  torque  for  substitution 
into  the  Euler  equations.  These  components  of  the  moment  vector  take  the  form 

LnFy(FF>y  +  LnFz<FF)z  +  LnRy(FR>y  +  LnRz(FR)z  (n  =  l-2'3) 
where  lnpy,  Lnpz,  LnRy,  and  LnRz  are  functions  of  the  known  parameters  of  the 
problem  only.  For  the  sake  of  brevity  in  writing  the  equations  of  constraint.- 
the  coefficients  are  given  here  as: 

llFy  "  em  cos6lx[sin  pjj  -  m  s(r)  cos  pjj] 

-  M's«»(r)  sin6lxcos  p^  +  (a-4'sr)  sindlx  sin  ^  ,  (49a) 

l1Fz  *  -8m  cos5lx[cos  p^  +  m's(r)  sin  pjjj] 

+  (a-p'sr)  sin6lx  cos  p^  +  p’sa(jj)  sin$ix  sin  ~  ,  (49b) 
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(49c) 


I  .  .  .  ns  ,  ,71  7TS 

LlRy  *  ®m  cosfilx  [sin  ^  cos  Fn] 

+  Sin6jx  cos  j~  -  (b+4'Br)  sin5lx  sin  ~  , 

i  .  r  ns  ,  ,7t.  ns. 

l1Rz  ■  ~em  ccs6lx[cos  ^  +  m'b^)  sin  pjj] 

-(b+M’Br)  sinfiixcos  ^  -  u'Bt>(j})  sin6lx  sin  ^  , 

L2Fy  ■  em  cos52x[sin  Ffi  -  M'S(J)  cos  pfi] 

+  M'sa(jJ)(::os5lx  cos<52x  cos  --5  +  sin  ~  cosa3x)/sin6lx 

(a-M'sr)(cos6ix  cos62x  sin  pfi  “  cos  ^  cos53X)/sin6lx  , 

l2Fz  *  ~em  cos62x[cos  pj-  +  n's(fj)  Sl'n  Fiji 

(a-ji'sr)  (cosSiy  cosfi2x  cos  ~  ♦  sin  ~  cos6iv)/sin6i v 

4'sa(jJ)(cos6lx  C0Sfi2x  sin  "  cos  ^  cos53x)/sin6lx  , 

L2Ry  x  *m  cos62x[sin  pjj  .  M'B(?)  cos  jjjj] 

~  ^’Bb(^) (costtjx  cosfi2x  cos  p^  +  sin  p--  cos$3x)/sin6ix 

.  JTS  7 T«5 

lD+M’BrJ  (cos6lx  cosa2x  sin  FR  “*  cos  FR  cos53x)/sin5ix  r 
l2Rz  *  -cm  cosfi2x[cos  pj  +  m'b(J)  sin 
(b+M'Br);cos5lx  cosfi2x  cos  pfc  +  sin  FR  cos53x)/sin6lx 
^’b^r)  (cos5jx  cosfi2x  sin  pjq  -  cos  pjq  cos53x)/sin6lx  , 


(49d) 


(49e) 


(49f) 


(499) 


(49h) 


THE  EQUATIONS  OF  CONSTRAINT 

We  are  now  in  a  position  to  write  the  six  equations  of  constraint.  This 
will  done  here  for  the  special  cases  of  a  perfectly  straight  tube  (ic*0),  and 
a  perfectly  balanced  and  aligned  projectile  (em  »  63X  *  0).  It  is  not  necessary 
to  make*  either  one  ov  these  simplifications  in  order  to  obtain  the  equations, 
but  by  making  one  or  the  other,  we  avoid  the  very  tediously  complex  equations  of 


the  general  case.  Furthermore,  this  arrangement  permits  the  study  of  the  iso¬ 
lated  effects. 

In  both  cases  it  will  be  assumed  that  the  only  motion  of  the  tube  is  that 
of  free  recoil  in  the  axial  direction.  Observing  that  (Qg)x  =  v*.  and  that  for 
free  recoil  the  only  forces  acting  on  the  tube  are  the  gas  pressure  force  and 
the  mutually  opposite  tube  to  projectile  forces,  we  have 

1  r r*t  ,"x  Bn,  *  ~  MBN.  .  _  .  _  _  . , 

vt  *  •  (r  ~  +  *  bfr  +  m  sff  -  pcai  <5° 

For  the  straight  tube  case  (k=0)  where  (Qg)y  =  (Qb)z  =  (“b'/Biy  = 

^ 

(“b'/BJz  *  (“b'/BJy  *  (Wb'/B)z  B  v< 

RB  *  SB  +  Pb  *  C(Qb)x  +  s]’ 

.its  ns  ,ns.a  ?rs 
“  mtrN  sin  rN  *  (rT5)  cos  rN  ]j 


Also, 


*  U  S 

rrs 

rN 

,  fiy  =  0 

,  flz  »  o  , 

. 

rrs 

*  w  = 

rN 

f  fly  =  0 

N 

IB 

o 

then: 

cos6^M^ 

The  moment  equations  are  then: 


+  LlFy(FF)y  +  l1Fz(ff)z  +  LlRy(FR)y  *  llRz(FR)z 

=  cos61x(^)It  *  cos62x  cos63x(p2)2(iT-iA)  , 
cosa2*M«  *  H(^g)[Pb*  *  -  Stm-— 

+  L2Fy(FF)y  +  L2Fz(Fp)z  +  L2Ry(FR)y  +  l2Rz(fr)z 
*  cosfi2x(^)Ir  -  cos63x  cos51x  (S)‘(IA.IT)  , 

and 

c°«3*mx  -  *  <£  -  -  ej 

+  L3Fy(FF)y  ♦  L3Fz(FF)z  +  L3Ry(FR)y  +  L3Rz(FR)z 
*  cos63x(^)lA  ■ 

For  the  balanced  and  aligned  projectile  case  (em  =  o3x  =  o)  where 

■***■*  -*  4  4 

(QB)y  ■  (Qb)z  *  («b'/B)y  =  (<*>bVB>z  -  (Wb’/B)y  =  (“b’/B)z  =  0  » 

RB  *  Sg  +  p0 

*  HQb)x  +  s]i  +  [k(s)*Ny]j  +  [k(s)*Ny]k  . 


The  force  equations  are  therefore 

<Pb  ♦  ?  Pel*  *  (1  ♦  ^)[<^  -  B»)(,t~MBl1) 
»t  c  mt'ur  N  ;Vb n  -  N> 

-  B  -  M'bfR  "  M’sff3  “  WpS  , 


(53a) 


(53b) 


(53c) 
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and 


(55a) 


(FF>y  ~  M*s(sHPf)2  *  <FR>y  *  M'etjjJWz 
*  n»pMs)*Ny  , 

<FF)z  +  M's^Wy  ♦  <FR>z  +  R>y 

«  IHpK(S)  *NZ  . 

Also,  53X  *  0  »>  S^x  *  $2x  *  90*  so  sinfiix  =  sin62x  *  cos63X  =  1,  and 

cos^ix  =  cosS2x  *  sin^3x  “  0* 

From  the  above  and  the  formulas  relating  the  direction  cosines.  Me  get 

„  ns  _  .ns 

cosfijy  a  COS  ~  ,  COSfil2  =  Sin  —  , 

.  .ns  „  ns 

COS62y  15  ~  Sin  *  Cos62Z  *  COS  PS  ' 

COSSjy  a  0  ,  COSfi3y  a  0  . 

The  moment  equations  therefore  become 

-♦  -4 

LlFy(FF)z  "■  LlFz(FF)z  ♦  LlRy(FRiy  ♦  LlRz*FRiz 

*  UKS  +  (^)a](Ny  sin  pfj  -  Nz  cos  p-j) 

-  KT(s)*(Ny  cos  ps  +  Nz  sin  ps> 

•.its.,.,  ns  „  .  ns ,  i  _ 

+  Ks(~)(Ny  cos  -  Nz  sin  p^) } It 

’/ns\/».  ns  ..  .  ns. . 

-  Ks(pjj)  (Ny  cos  ♦  Nz  sin  p^HIt-Ia)  r 

L2Fy(FF)y  +  L2Fz(FF)z  +  L2Ry(FR)y  +  L2Rz<FR>z 

*  {(kS  +  55  (s)*](Ny  cos  p^  +  Nz  sin  pjq) 


+  ict(s) * (Ny  sin  pjjj  -  Nz  cos  ~) 

-  KS(-)(Ny  sin  ^  +  Nz  cos  -)}IT 
“  *s(p|)(Ny  sin  PR  _  Nz  cos  PrHiT"ia)  • 


(56a) 


(56b) 
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and 

Mx  +  *-3Fy (Fp)y  +  L3FZ(Ff)z  +  l-3Ry(FR)y  +  L3Rz(FR)2  =  —  Ia  •  (56c) 

A  brief  note  and  suggestion  concerning  the  solution  of  the  equations  would 
seem  to  be  in  order  here.  In  the  first  force  equations,  there  appear  the  terms 
Fr  and  Fp  which  as  nonlinear  functions 

FR  -  ^fRly  *  <FR>i  . 

and 

FF  =  ’'Wy  +  (Pp)z 

of  the  unknowns  would  seem  to  present  a  difficulty.  This  difficulty  should  be 
fairly  easy  to  circumvent  by  first  substituting  estimates  for  FR  and  Fp  and 
solving  the  equations  for  S.  Estimates  of  FR  and  Fp  can  be  based  on  their  imme¬ 
diately  previous  up-bore  calculation  (FR  =  Fp  =  0  st  s  =  0)  and,  as  inspection 
of  the  equations  show,  are  not  critical.  Then,  after  integrating  to  get  s 
and  substitution  of  its  value  and  that  of  s  where  required  in  the  remaining 
equations,  the  simultaneous  solution  for  (FF)y,  (Fp)z,  (FR)y,  (FR)Z,  and  Mx  is 
obtained.  FR  and  Fp  are  then  calculated  and  compared  to  the  original  estimate. 
If  too  different  from  the  estimate,  the  calculated  values  can  be  used  and  the 
above  repeated.  Following  this  solution  outline,  the  equations  hava  been 
programmed  for  computer  solution  for  a  sp  ctrum  of  gun  and  ammunition  com¬ 
binations,  as  will  be  reported  in  a  subsequent  report. 

Also,  it  should  ba  noted,  in  arriving  at  these  demonstration  equations,  the 
terms  that  couple  the  dynamic  response  of  the  tube  and  the  constraint  forces 
were  dismissed  by  assumption.  In  reality,  of  course,  the  force  of  the  projec¬ 
tile  bearing  against  the  tube  causes  motion  and  flexure  of  the  tube  which  in 
turn  have  an  effect  on  the  forces  which  are  their  cause.  For  some  studies  this 
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nay  not  be  inportant.  But  for  others,  it  would  be  a  requirement  to  solve  the 
equations  of  constraint  and  the  tube  dynamics  equations  simultaneously.  This 
suggests  the  utility  the  equations  of  constraint  may  have  in  tube  dynamics 
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APPENDIX 


DERIVATION  OF  THE  DIRECTION  COSINE  RELATIONSHIP 
Since  we  shall  be  concerned  here  only  with  axially  symmetric  projectiles, 
it  will  be  permissible  to  arbitrarily  establish  the  initial  position  of  one  of 
the  transverse  axes  of  inertia  with  respect  to  the  b  coordinate  system.  Thus, 
the  number  of  rotational  degrees  of  freedom  of  the  body-fixed  coordinate  system 
with  respect  to  the  b  coordinate  system  is  reduced  from  three  to  two.  Because 
of  the  simplification  in  analysis  that  is  affected,  we  shall  select  the  1-axis 
to  lie  initially  in  the  x-y  plane  as  shown  in  Figure  A-l. 


Figure  A-l.  Initia1  position  of  the  body-fixed  coordinate  system  in 
relation  to  the  moving  coordinate  system. 

The  direction  cosines  of  the  initial  position  of  the  1-axis  are  related 

according  to 


COS*6*ix  +  COS*5*iy  +  003*5*12  *  1  .  (A-l) 

Now  since  6*ix  *  fiix  *  constant,  and  6*iz  *  90*,  we  have 

cos*6*iy  =  l  -  cos*5ix  *  sin*5ix  .  (A-2) 

Figure  A-2  shows  an  axis  parallel  to  the  1-axis  and  an  axis  originally 
parallel  to  the  y-axis,  but  fixed  in  the  body  and  therefore  rotated  through  the 
angle  6,  the  angle  of  rotation  of  the  projectile. 
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Figure  A-2.  Auxiliary  axes  fixed  in  the  body. 

In  the  figure  the  plane  defined  by  triangle  AOC  is  perpendicular  to  the 
plane  defined  by  triangle  BOC,  and  the  plane  defined  by  triangle  ABC  is  perpen¬ 
dicular  to  line  0y<.  Therefore,  the  angles  ACB,  ACO,  and  BCO  are  right  angles 
and  if  we  set  GC  =  1,  then  CA  *  tan©,  BC  *  tan6ejy,  OA  *  1/cos©,  and  OB  = 
l/cos5*jy.  From  the  law  of  cosines 

cos«ly ,  <”1:  -  <“>■ .  (A-3 

2 (OB) (OA) 


Now  since 


and 


tan1© 


tan*5*^y 


1  -  cos*© 

COS  *  (? 


I  -  cos16*1y 

COS*6*Jy 


we  get 


cos© 


ly 


COS*©iy  cos©  1  -  1 


cos 


+  COS*5*ly  1 

- jr - 1  *  ~ 


1  •+  cos*© 


ly 


COS*  0 


] 


■  cosd'jy  cos©  =  t  sinfijx  cos©  . 

The  positive  sign  is  appropriate  here  and  so  for  constant  twist  rifling 

cosfijy  n  sin6jx  cos  j-jj  . 


(A-4) 


(A-5) 
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Also 


7TS 


cos*6jz  »  1  -  cos*6jx  -  sin*6jx  cos*  pjj 


so  that 


cos5 


12 


.  .  .ns 
sinSlx  sin  ^ 


(A-6) 


(A- 7  ] 


The  four  remaining  direction  cosine  variables  can  now  be  determined  with  two 
equations  of  the  normalization  condition  and  two  equations  of  the  orthogonality 
condition: 


COS*52x  +  COS*62y  +  COS*52Z  = 

1  , 

(A-0a) 

cos*63x  +  cos*63y  +  cos*63z  = 

1  , 

(A-8b) 

cos5jx  cos62x  +  Cos ly  COs62y  +  cos5iz 

COSd2z  =  0  , 

(A-9a) 

COSfiix  COSd3x  +  COSfijy  COS63y  +  C0sfii2 

cosfi3z  *  0  . 

(A-9b) 

Solving  for  cosd2z  in  Eq.  (A-9a)  gives 

coso2z  ■ 


COSdly  COS$9y  +  COSdlw  COSfi^w 


COS6iz 

Substitution  into  Eq.  (A~8a)  and  rearrangement  then  gives 

cos*5ly  C0S*52y  +  cos*$2y  cos*6i2 
+  2cos6ix  cosd2x  cosily  cos62y 
♦  COS*$2X  COS*djz  +  COS*5jx  COS*52x  -  COS*$iz  =  0  . 


(A-10) 


( A— 11) 


T  f  wa  1  of  a  s  rrte  t  A . 


a  * . 


—  c 


fiy  '  V1Z '  w  *"  a  mu  u  —  l,u:»-c/2x 

cos*5jz  +  cos*^*  +  cos*62x  "  cos*5jz,  then  by  the  quadratic  formula 

■b  ±  /b*-4ac 


COS52y  * 


2a 


cosfi 


2y 


-COS6^x  COS62X  cosily  ±  cosg^/cos*  6jz+COS*  i?2y-COS*62x 


cos*5ly  +  cos *6 
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(A-12J 


Now  since 


C°S*«1X  ♦  cosa*ly  +  cos5*jz  =  1  o  cos*5lx  +  cos*«2x  +  cos*<53} 


I 
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we  have 

cosf5iz  +  cos*$iy  -  cos*62x  ~  cos*63X  .  (A-13) 

The  expression  under  the  radical  in  Eq.  (A-12)  is  thus  simplified  so  we  now  have 


cosS 


2y 


-  cos6jx  cos62x  cos6|y  i  cos6iz  cosd3X 


cos*6ly  +  cos?aiz 

Substitution  into  the  expression  for  cos62Z  and  rearrangement  now  gives 


(A-14) 


-  cos5jx  cos62x  cosalz  ±  coS(5iycos53x 

COS622  s  — — ■■  —  . . * .  (A-15) 

cos*6^y  +  cos*6lz 

The  t  signs  in  these  expressions  for  the  direction  cosines  can  be  resolved 
by  expressing  mathematically  the  fact  that  the  1,  2,  3,  coordinate  system  is  a 
right-handed  cartesian  system.  That  is,  if  the  1  axis  is  crossed  into  the  2- 
axis  the  result  is  the  3-axis. 

-•  -4  ^ 

(cos6jxi  +  cos6lyj  +  cos5lzk)  x  (cos^x-*  ♦  cos62yj  +  cos62zk) 

^  -f  -♦ 

=  cos6lx  cos62yk  -  cos5jx  cos«2zJ  ~  cos6jy  cosfi2xk 
+  cosfijy  cosd2zi  +  cos6jz  cos$2xJ  “  cos6jz  cos^yi 

-4  ^ 

*  cosfi3xi  +  cos«3yj  +  cos63zk  .  (A-16) 

Equating  the  coefficients  of  the  same  unit  vectors  of  both  sides  of  the  equation 
gives  for  the  i  unit  vector 

CvS5jy  COSu2Z  “  COsoj2  COSo2y  =  COS63x  .  ( A- i 7 ) 

Substitution  of  Eqs.  (A-14)  and  (A-15)  and  rearrangement  gives  the  following: 
-cosfijy  cos$ix  cosfi2x  cosfij-  ±  COS2fijy  cos$3x 
+COS6jz  COSfijx  COSd2x  cosily  T  cos*<5jz  cos63x 

■  cos63x(cos*6iy  +  cos*6i7)  .  (A-18) 

For  equality  to  exist  in  the  above  equation,  it  is  clear  that  the  signs 
must  be  as  follows: 


-cos5jx  cos52x  cos6jy  -  cos$iz  cos63x 
cos52y  “  cos*fiiy  +  cos*6jz 


(A-19) 
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and 


-C0S6ix  CCS02X  cos6iz  +  COS6lv  COS53x 

cos 5o?  ■  —  ~  ~  ..  —  < -  . 

££  cos*fiiv  +  cos*6iz 


>ly  +  cos^fij^ 
In  the  same  way  as  above,  it  can  be  shown  that 


-cos$ix  cos63X  COS^ly  +  COSOI2 

COs53y  *  cos*6iy  +  cos46iz 


cosfi2x 


and 


cos63Z 


-cos6ix  cosfi3X  cos5iz  -  cosfijy  cos62x 

^  COS'^iy  +  C0S*6iz 


Now  since 


and 


cos6ly  *  sin«lx  cos  p^ 


cos5lz  *  sin$lx  sin  -r. 


for  constant  pitch  rifling  we  also  have 

cos*6ly  +  cos*$j2  »  sina6ix(cos8  p^  +  sin*  pj^) 


so  that  finally  we  have  the  direction  cosine  relationships  shown  as  Eqs. 
through  (26f). 


(A-19b) 

(A-20) 


(A-21 ) 
(26a) 
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rB  =  SB  +  PB 
Texts  on  dynamics  (ref  10)  g->ve  the  formula 

■  >  •  •  • 

-*—*—*•  — *■  — f  -¥ 

pB  =  pb  +  o>b/B  x  pb  +  2wb/B  x  pb  +  wb/B  x  <ob/B  x  pb  (6) 

for  rotating  coordinate  systems  where  the  subscript  b/B  indicates  rotation  of 
the  b  coordinate  system  with  respect  to  the  B  coordinate  system.  If  em  is  the 
magnitude  of  the  eccentricity  of  the  projectile  c.g.  with  respect  to  the  spin 


axis,  then  pb,  pb,  and  pb  can  be  expressed  as 

Pb  =  em  (c°s6 j  +  sin0k)  , 


and 


pb  *  em8  (-sinej  +  cosGk.)  , 


pb  =  em[-(G  sin0  +  0 2  cos0)j  +  (0  cosO  -  02  sinG)k] 
where  Q  is  as  shown  in  Figure  2. 

V 


(7) 

(7a) 

(7b) 


Figure  2.  Position  vector  for  the  projectile  c.g.  in  the  moving 
coord  i  n3to  systi sm . 

Since  the  spin  axis  is  taken  in  general  to  be  the  tangent  to  a  curve  moving 
in  space,  it  will  be  convenient  to  refer  to  the  Frenet-Serret  Formulas  of  dif¬ 
ferential  geometry.  These  are 


■p  =  -rN  ,  =  kN  ,  and  ;p  =  tB  -  kT 

ds  ds  ds 


(8a,b,c) 


o 


where 


—  -»  -»  -*■  (-IR  1 

B  =  T  x  N  ,  T  =  "  ,  and  k.  ^  . 

ds  p 

The  velocity  of  a  point  moving  along  a  curve  in  space  is  expressed  as 


v(t)  = 


dR  _  ds  dR 
dt  "  dt  ds 


Making  use  of  the  Frenet-Serret  formulas  allows  the  velocity  to  be  written  as 


and  the  acceleration  to  be  written  as 


H2c  -*■  He  2  -* 

a(t)  *  h  T  *  <a;>  KN  • 


Now  since  Sq  is  the  sum  of  the  acceleration  of  the  origin  of  the  b  coor¬ 
dinate  system  relative  to  its  path  and  the  acceleration  of  the  point  on  the  path 
coincident  to  the  b  coordinate  system  origin,  we  can  write 

sB  *  QB  +  Hi  T  +  (jjjf)*icN  .  (12) 

Also  note  that  %/b  "is  the  sum  of  the  angular  velocities  due  to  rotation  of  the 
spin  axis  at  the  point  coincident  with  the  origin  of  the  b  coordinate  systam  and 
motion  of  the  coordinate  system  origin  relative  to  its  path.  Therefore,  if  we 
designate  by  the  subscript  p:,  a  coordinate  system  which  is  coincident  to  the  b 
coordinate  system  but  fixed  to  the  spin  axis,  we  can  write 

— ►  -*•  -4 

Wb/B  =  Wb'/B  +  wb/b'  •  U3) 

Recalling  that  wp/g  has  no  component  in  the  direction  of  the  spin  axis,  it 


follows  that 


T,  ds;  dT  ds!  ~ 

ub/b'  =  dt  T  x  di  =  dt  T  x  KN 


6 


so  that 


“b/B  =  «b*/B  +  al  T  X  KN 


ds  p 
“b'/B  +  3t  kB 


Differentiating  with  respect  to  time  gives 


Wfo/B  =  Wb'/B  +  5i;  [5^  KB1 


which,  after  some  manipulation,  results  in 


r  d2s  die  ,ds  2  r  ^  .ds 
"b/B  =  wb ’ /B  +  U  3^  +  0i  (dt}  38  "  KT(dt)N 


COORDINATE  TRANSFORMATION 

First  of  all,  by  observing  that  T  s  i,  and  N  =  Nyj  +  Nzk,  it  follows  that 
-*■  -♦ 

B  =  T ‘N  =  Nyk  -  Nzj.  Also, 


Therefore, 


Qb  =  (QB)xi  +  (QB)yj  +  (QB)zk  » 

(18a) 

-f  -♦  -»  -♦  — *■ 

Wb'/B  *  (wb'/B)yJ  +  (“b'/Bizk-  > 

(18b) 

•  •  % 

— * 

"b'/B  "  («D'/B)yj  +  (Wb'/B)zk  • 

(18c) 

SB  =  [<Qb)x  +  Hill  *  t(QB)y  +  (a?)J,cNy]j 
+  l(9b)z  +  (a?)2KNz]k  . 


ds  ..  .T  ds 

“b/B  "  n«b!/B>y  -  3t  kn2jj  +  L^b'/B )z  +  dt  KNyJK 


(19a) 

(19b) 


r , .  .  d2s  die  .ds,2!..  ,ds *  »T 

«b/B  =  {(«b'/B>y  "  t*  at2  +  di  (dt>  3Nz  "  ^dt3  Ny»J 

,  ,  .  d2s  die  .ds*,*.  ,ds,*u  ip 

+  {(«b'/B)z  +  at2  +  di  (dt3  ]NV  ~  KT<dt)  Nz’k  ■ 


(19c) 


The  total  angular  velocity  of  the  projectile  is  the  vector  sum  of  the  pro¬ 
jectile  spin  velocity  and  the  angular  velocity  of  the  b  coordinate  system,  or 

n  =  W  +  Ub/B  •  (2°) 
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In  the  b  coordinate  system,  the  components  of  are  then 

=  ^  < 

-  ds 

Oy  =  (Wb'/B)y  +  KNz  ■  and 
=  (Wb'/B)z  +  5t  K?'y  ' 

Similarly,  the  angular  acceleration  and  its  c<  mponents  are  given  by 


(21a) 

(21b) 


v*- 


lc) 


and 


1  =  w 

+  wb 

3  • 

(22) 

«x 

=  a) 

r 

(23a) 

d2s 
dtz  + 

d< 

ds 

^dt 

^z 

ds  2 

-  KT<3?>  Ny  ' 

(23b) 

d2s 
dt2  + 

dx 

ds 

(~) 

'dt' 

Z]Ny 

ds  2 

-  KT(^)  Nz  . 

(23c) 

It  remains  now  to  express  the  projectile  angular  velocity  and  acceleration 
in  terms  of  their  components  in  the  body-fixed  coordinate  system.  For  the  angu¬ 
lar  velocity 

=  cosoXxQx  +  cosfijySly  +  cos6Xz(2z  ,  (24a) 

fl2  =  C0S52x^x  +  cos62yfly  +  cos62z7lz  ,  (24b) 

and 

83  =  cos63Xilx  +  cos63yQy  +  cos63Zilz  .  (24c) 

Remembering  that  <5Xx,  d2x.  and  <53X  are  assumed  to  remain  constant,  we  get  by 
differentiation 

fij  =  COScSlxfix  +  (COS6ly)Qy  +  COSdjyQy 

d  , 

+  ^  (cos5lz)nz  +  cos6Xzflz  ,  (25a) 

*  •  d  • 

n2  =  COS62xtlx  +  (C0S<52y)£2y  +  COSfi2y^y 

+  ~  (cos62z)?iz  +  cos52zQz  ,  (25b) 
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